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Abstract 


This  project  has  demonstrated  that  one  class  of  magnetic  pulsations  known 
as  stormtime  Pc  5  waves  is  correlated  with  substorm  onsets.  Siorrntime  IV  5 
waves  observed  by  geostationary  satellites  in  the  afternoon  sector  is  character¬ 
ized  by  oscillations  of  magnetic  field  with  a  period  from  2  to  !0  minutes,  easily 
detected  by  magnetometers  on  communication  or  weather  satellites.  The  es¬ 
timated  substorm  unset  times  are  found  to  be  within  20  minutes  of  the  actual 
substorm  onset  times.  Geosynchronous  satellites  in  the  afternoon  sector  would 
detect  these  low  frequency  wave  events  about  2-4  hours  after  a  substorm  on¬ 
set  occurring  at  local  midnight.  The  delay  time  depends  on  the  propagation 
velocity,  which  varies  from  a  few  km/'s  up  to  50  km/s.  The  disturbed  region  of 
a  stormtime  Pc  5  event  has  a  longitudinal  extent  varying  between  30  and  00 
degrees.  The  study  shows  that  stormtime  Pc  5  waves  have  a  wave  amplitude 
confined  with  about  10°  from  the  magnetic  equator.  The  propagation  velocity 
is  found  to  increase  with  wave  frequency  and  with  the  magnetic  field  inclina¬ 
tion  angle.  Comparison  of  the  statistical  properties  of  stormtime  Pc  5  waves 
with  theoretical  calculations  of  propagation  velocity  suggests  that  the  propa¬ 
gation  velocity  of  storrntime  Pc  5  waves  agrees  better  with  the  perpendicular 
group  velocity  of  drift  mirror  mode.  The  propagation  velocity  of  stormtime  Pc 
5  waves  appears  to  be  mainly  determined  by  wave  parallel  wavelength,  which 
is  in  turn  determined  by  the  inclination  angle  or  the  magnetic  field  topology 
The  obtained  results  about  the  propagation  properties  of  magnetic  pulsations 
during  storm  times  is  important  for  the  satellite  operation  since  it  can  be  used 
to  predict  the  plasma  environment  a  synchronous  satellite  might  encounter. 
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I  Introduction 

This  report  summarizes  the  research  findings  obtained  during  the  support  of 
the  Air  Force  Office  ot  Scientific  Research  ( Al'  OSR  )  grant  F  19G20-S9-(  ’-000S. 
One  objective  of  the  grant  is  to  investigate  the  feasibility  of  using  dual  satel¬ 
lite  observations  of  low  frequency  magnetic  pulsations  to  remotely  diagnose 
substorm  onsets.  Another  objective  of  the  study  is  to  understand  tin*  propa¬ 
gation  velocity  of  stormtime  Pc  5  waves,  including  its  dependence  on  plasma 
variables  and  its  correlation  with  physical  parameters. 

Synchronous  satellites  stationary  over  North  America  generally  detected 
low  frequency  magnetic  pulsation  about  2-4  hours  after  substorm  onsets  oc¬ 
curring  near  local  midnight  over  the  Russian  stations.  In  this  project  we  in¬ 
vestigated  the  feasibility  of  remotely  sensing  auroral  substorm  activities  near 
midnight  by  using  magnetic  field  measurements  from  synchronous  satellites  on 
the  dayside.  A  synchronous  satellite  observing  magnetic  activities  on  the  day- 
side  's  difficult  to  measure  remotely  auroral  ubstorm  activities  occurring  on 
the  other  side  of  Earth.  The  knowledge  of  substorm  activities  is  important  for 
the  satellite  operation  since  it  can  he  used  to  predict  the  plasma  environment 
a  synchronous  satellite  might  encounter.  During  auroral  substorm  activities, 
energetic  plasma  has  been  injected  into  the  nightside  magnetosphere,  and  pre¬ 
cipitated  into  the  iono  phere  along  magnetic  field  lines,  producing  auroras, 
ionospheric  disturbance,  and  electromagnetic  radiation  over  a  wide  frequency 
range.  As  a  result,  spacecraft  are  sometimes  charged  to  a  higher  potential, 
affecting  the  operation. 

GOES  synchronous  sa  ellite  magnetic  field  data  since  1979  have  been 
surveyed  and  statistically  correlated  with  magnetograms  from  Russian  ground 
stations  in  this  project.  It  is  found  that  one  class  of  magnetic  pulsations 
known  as  stormtime  Pc  5  waves  is  correlated  with  substorm  onsets.  This 
class  of  magnetic  pulsations  is  characterized  by  oscillations  of  magnetic  field 
with  a  period  from  2  to  10  minutes,  easily  detected  by  magnetometers  on 
communication  or  weather  satellites.  Because  these  waves  propagate  from 
midnight  to  dayside  after  substorm  onsets  with  high  speed,  they  can  be  used  to 
remotely  diagnose  substorm  activities  when  geosynchronous  satellites  observe 
them  on  the  dayside. 

Figure  1  illustrates  the  scnematic  of  determining  the  propagation  velocity 
of  magnetic  pulsations  from  GOES  magnetic  field  measurements.  Two  GOES 
satellites  usually  separated  by  two  hours  will  detect  the  onset  of  a  low  fre¬ 
quency  magnetic  pulsation  event  within  a  period  ofl0-30  minutes.  As  shown 
in  the  schematic,  the  onset  time  at  each  GOES  satellite  is  determined  by  the 
peak  of  the  first  oscillation.  The  propagation  velocity  is  then  calculated  from 
the  satellite  separation  distance  D  and  the  time  difference  A t  between  the 
onset  timesobserved  by  the  two  GOES  satellites.  Assuming  a  substorm  on¬ 
set  near  local  midnight,  the  substorm  onset  time  is  estimated  from  the  wave 
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propagation  velocity  and  the  distance  >  horn  local  midnight  to  the  satellite 
( Figure  1  ). 

I'lie  results  obtained  in  this  study  suggest  that  stormnme  Pr  5  waves 
are  correlated  with  subsiorin  onsets.  The  estimated  substorrn  onset  times  are 
found  to  be  within  20  minut*-  of  the  actual  substorm  onset  times.  Geosyn¬ 
chronous  satellites  in  the  afternoon  sector  would  detect  these  low  frequency 
wave  events  about  2-1  hours  after  a  substorm  onset  occurring  at  local  mid¬ 
night.  The  delay  time  depends  on  the  propagation  velocity,  which  varies  from 
a  few  km/s  up  to  50  km/s.  The  propagation  velocity  is  found  to  increase  with 
wave  frequency,  which  is  found  to  be  correlated  with  the  magnetic  field  inch 
nation  angle.  Comparison  of  the  statistical  properties  of  stonntime  Pc  5  waves 
with  theoretical  calculations  of  propagation  velocity  supports  that  stormtime 
Pco  waves  propagate  as  the  drift  mirror  mode.  The  propagation  velocity  of 
stormtime  Pc  5  waves  appears  to  be  mainly  determined  by  wave  parallel  wave¬ 
length.  which  is  in  turn  determined  by  the  inclination  angle  or  the  magnetic 
field  topology. 

Below  we  briefly  summarize  the  results,  and  describe  in  detailed  the  anal¬ 
ysis  in  the  following  sections. 

A  Remote  Diagnostic  of  Substorm  Onsets 

We  conducted  data  analyses  of  GOES  dual  satellite  magnetic  field  data  for  the 
years  of  1979.  1983.  and  1986.  We  investigated  the  correlation  of  stormtime 
Pc  o  waves  with  substorm  onsets.  From  the  data  set.  we  selected  stormtime 
Pc  5  events  that  ground  station  magnetometers  were  available.  Because  most 
correlated  ground  station  magnetometer  data  were  collected  in  Russia,  we 
were  able  to  conducted  the  correlation  study  only  for  the  data  set  of  1979.  We 
described  the  results  in  Section  II. 

This  study  used  simultaneous  observations  of  stormtime  Pc  5  events  by 
the  GOES  2  and  GOES  3  geostationary  satellites  to  study  the  correlation  of  Pc 
5  waves  with  substonn  onsets.  Eighteen  Pc  5  events  occurring  from  March  to 
December  1979  were  first  surveyed.  After  excluding  events  with  highly  vary¬ 
ing  activity  or  scarce  station  coverage,  only  six  events  with  clear  substorm 
features  2-4  hours  before  the  GOES  observation  of  the  Pc  5  events  were  then 
analyzed.  From  the  wave  propagation  speed  and  the  distance  from  local  mid¬ 
night  to  the  satellite  position,  the  substorm  onset  time  was  estimated.  For 
the  six  events  examined,  ground  magnetograms  near  local  midnight  indicate 
that  a  substorm  onset  occurred  within  20  minutes  of  the  estimated  substorm 
onset  times.  This  result  suggests  that  the  occurrence  of  stonntime  Pc  5  waves 
is  probably  correlated  with  substorm  onsets.  The  results  were  published  in 
Journal  of  Geophysics  Research  (Pangiaet  al..  1990). 
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Figure  1:  Schematic  of  determining  propagation  velocity  of  low  frequency 
magnetic  pulsations  by  using  dual  GOES  satellites. 
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B  Statistical  Survey  of  Propagation  Velocity  of  Storm- 
time  Pc  5  Waves 

We  conducted  a  statistical  study  of  propagation  velocity  of  storm  time  Pc  5 
waves  observed  by  dual  GOES  satellites  during  the  years  of  1979.  IDS.'i.  and 
1986.  We  tabulated  the  properties  of  these  events  and  deduced  the  propa¬ 
gation  velocity  from  the  conjunction  observations,  the  GOES  satellite  data 
was  surveyed  to  deduce  the  seasonal  variation  for  the  occurrence  of  Pc  b  wave 
events  and  the  distribution  of  propagation  velocity  over  frequency  and  local 
time.  We  also  conducted  a  correlation  study  to  find  the  relationship  between 
the  propagation  velocity  and  wave  properties.  It  is  found  that  the  propaga¬ 
tion  velocity  increases  with  frequency,  which  increases  with  the  magnetic  held 
inclination  angle.  The  variation  of  propagation  velocity  on  parameters  was 
also  used  to  determine  the  wave  mode  responsible  for  storm  time  Pc  b  waves. 
These  results  were  presented  in  the  1992  Spring  AGE  meeting  (Lin  et  ah. 
1992).  These  results  are  described  in  Section  III. 

C  Theoretical  Study  of  Low  Frequency  Wave  Propa¬ 
gation  Velocity 

In  addition  to  data  analyses,  we  modeled  the  wave  propagation  velocity  to 
better  understand  the  statistical  properties  of  stormtime  Pc  b  waves  and  the 
correlation  bewteen  magnetic  pulsations  and  substorm  onset.  We  solved  the 
wave  dispersion  equation  and  calculated  the  wave  group  velocity  for  several 
wave  modes  that  have  been  suggested  to  be  responsible  for  exciting  stormtime 
Pc  5  waves.  The  numerical  results  are  given  in  Section  IV. 

D  Eigenmode  Analysis  of  Low  Frequency  Waves  in  a 
Dipole  Magnetic  Field 

Since  the  Earth's  dipole  magnetic  field  could  limit  wave  amplitude  structure 
along  field  lines  and  thus  affect  the  synchronous  satellite  observations  of  mag¬ 
netic  pulsations,  we  also  performed  an  eigenmode  analysis  of  low  frequency 
waves  in  a  dipole  magnetic  field.  The  study  shows  that  stormtime  Pc  5  waves 
have  amplitude  confined  within  about  10°  from  the  magnetic  equator.  This 
result  explains  that  GOES  satellites  rarely  detected  stormtime  Pc  b  waves 
during  summer  months.  The  calculations  are  described  in  Section  V. 


•1 


II  Remote  Diagnostic  of  Substorm  Onsets 

Stonutinif-*  Pc  5  waves  have  long  been  known  to  be  associated  with  eeomat»- 
netic  activity  (Bartield  and  Coleman.  1970:  Barfield  et  ai..  1972:  Barfield  and 
McPherron.  D72).  Typically,  they  occur  during  the  main  phase  of  a  genuine 
lU'tic  storm  iti  the  afternoon  sector.  Barfield  and  McPherron  i  1972}  lompaied 
ATS  1  observations  of  Pc  5  events  to  ground  magnetograms  and  found  that 
seventeen  of  their  twenty  events  studied  were  closely  correlated  with  the  on 
set  of  a  substorm.  The  objective  of  this  study  is  to  further  investigate  this 
correlation  with  mult iple  satellite  observations  of  Pc  5  waves. 

Multiple  satellite  observations  olfer  additional  information  over  single 
satellite  observations  regarding  propagation  characteristics  oi  Pc  ."events  (Al¬ 
lan  et  ah.  1982:  Walker  et  ah.  1982:  lakahashi  et  ah.  1985;  l.m  and  Barfield. 
1985).  for  example,  a  study  of  global  t ompressiuna!  Pe  5  wave  observations 
by  Takahashi  et  al.  (198.5)  indicates  that  the  waves  propagate  westward  m 
the  afternoon  sector  with  speeds  in  the  range  of  !  il  kites.  In  a  statistical 
study  of  stormtime  Pc  5  waves.  Tin  and  Barfield  (1985)  found  that  the  wave 
propagation  speed  is  typically  less  than  30  km/s. 

For  this  study,  propagation  information  was  used  as  an  additional  test  of 
i ho  correlation  between  stormtime  Pc  5  waves  and  substorm  onsets.  By  using 
the  azimuthal  propagation  velocity  of  stormtime  P<  5  wave's,  we  first  estimated 
the  substorm  onset  times  and  then  correlated  the  estimated  onset  times  with 
the  actual  substorm  onset  times  determined  from  the  ground  magnetogram 
data. 

A  Method 

The  propagation  velocities  of  stormtime  Pc  5  waves  in  this  study  wen*  taken 
from  Lin  and  Barfield  ( 1985).  who  examined  30  Pc  5  events  simultaneously  ob¬ 
served  by  GOES  2  and  GOES  3  satellites  during  the  one  year  interval  March 
1979  to  February  1980.  During  the  study  interval,  the  two  .satellites  wen- 
located  within  approximately  2  hours  of  local  time  (approximately  30°  lon¬ 
gitude)  of  each  other.  Both  satellites  were  on  the  geographic  equator,  with 
GOES  2  at  104°  ±  1°  west  geographic  longitude  and  GOES  3  at.  135°  west 
geographic  longitude.  I  lie  onset  times  of  Pc  5  wave  events  at  each  satellite 
Tw  were  determined  to  be  at  the  peak  ot  the  first  magnetic  field  oscillation. 
The  propagation  velocity  c7  is  men  determined  from  the  satellite  separation 
distance  D  and  the  time  difference  between  Pc  5  wave  event  onset  times  G/\, 
as 

v.j  =  D/A7V-V,  U) 

where  V,  =  3.09  km/s  is  the  satellite  velocity  at  the  synchronous  orbit.  The 
error  in  deducing  the  wave  propagation  was  found  to  be  mainly  due  to  uncer 
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taint  ies  in  identifying  the  tirst  oscillation  since  flu*  wave  amplitude  delected 
i>Y  the  satellites  sometimes  mows  Ironi  <t  grading  depression  ot  the  magnetic 
field,  i  here*  the  <-nor  <>(  the  onset  time  would  he  at  most  one  wave  pejiod. 
1’sinit  tit:.-  lefinitiou.  we  found  that  the  error  in  estimating,  the  propagation 
veloci  ,s  generally  less  than  i  Vi . 

lo  estimate  suhstonn  onset  times  from  IV  j  propagation  speeds,  we  as- 
siime  that  (1/  plasma  injection  occurs  near  midnight  during  suhstonn  onset. 
(2)  the  injected  plasma  excites  stormtime  Pc  h  waves  observed  by  GOES  satel¬ 
lites  in  the  afternoon  sector,  and  (3)  the  Pc  5  waves  propagate  at  t  lie  average 
drift  speed  of  the  injected  plasma.  These  assumptions  are  based  on  the  satel¬ 
lite  observation  reported  by  l.anzerotti  el  al.  (197Y)  that  1  HI  ke\  protons  are 
associated  with  a  stormtime  Pc  wave.  Furthermore,  our  study  neglects  the 
radial  propagation  velocity,  which  may  sometimes  be  appreciable.  However, 
no  observation  of  the  radial  propagation  has  been  reported. 

Since  (.in  and  Barfield  (  1 98  d  j  had  already  surveyed  the  (JOBS  magnet  i< 
held  data  for  stormtime  Pc  d  wave  events  occurring  during  the  one  year  period, 
we  decided  to  survey  ground  magnetometer  data  during  these  wave  events 
Among  the  SO  wave  event"  in  I  he  dataset  of  l.iti  and  Hatfield  I  P.l'vdf.  we  were 
able  It)  obtain  ground  magnetometer  data  near  midnight  tor  Is  events,  We 
estimated  tin*  time  difference  between  the  substorm  onset  and  tin*  GOES  ob¬ 
servations  of  the  stormtime  Pc  d  waves  by  dividing  the  distance  from  GOES  2 
to  local  midnight,  by  the  wave  group  velocity.  Extrapolating  the  wave  events  to 
local  midnight  shows  that  magnetograms  from  Russian  ground  stations  were 
required  to  conduct  this  survey.  After  analyzing  these  ground  magnetograms, 
we  found  that  magnetic  activity  was  generally  present  for  all  the  events  dose 
to  the  predicted  substorm  onset  times.  However,  some  ground  magnetograms 
were  not  suitable  for  study:  we  excluded  seven  events  from  the  study  because 
of  highlv  varying  activity,  and  another  five  events  because  of  scarce  station 
coverage.  We  based  our  study  on  the  six  events  that  showed  clear  substouri 
onsets  (listed  in  Table  1  |. 

B  Data  Analysis  Results 

We  present  two  examples  ot  COES  and  ground  magnetic  field  data,  figure  2 
shows  the  three  components  of  the  GOES  magnetometer  data  for  a  two  hour 
period  starting  at  21:30  f  l  on  March  26.  PE 9  when  GOES  2  was  at.  about 
1  1:30  local  time.  The  orientation  of  the  components  are  that  lip  is  parallel  to 
the  Earth’s  rotation  axis.  He  is  in  the  radial  direction,  and  Hn  completes  the 
orthogonal  coordinate  system.  A  Pc  5  wave  event  with  about  a  seven  minute 
period  is  seen  in  the  lip  and  He  components  of  the  GOES  3  measurements 
from  22:16  to  about  22:10  FT.  The  GOES  2  satellite  detected  the  wave  before 
GOES  3.  indicating  a  westward  propagation  since  GOES  2  was  eastward  of 
GOES  3  bv  30°  in  longitude.  At  GOES  3.  Hie  onset  time  of  the  wave  event 
was  determined  to  be  22:16  FT  at  the  first  peak  in  the  wave  train  of  the 
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TABLE  1.  Comparison  of  Actual  and  Estimated  Substorm  Onset  Times 


t  ase 

Date 

7V 

(UT)  (MIT) 

i  km/sec) 

(UT) 

rJ 

t  UT) 

l  it  --  />- 
i  mm  1 

Te  - 

l  mu 

Ts 

) 

1 

79-03-25 

21:54 

14:56 

16.3  ±  2.0 

20:11  ±  15 

19:55 

119 

1 6  ± 

!  5 

*> 

79-03-28 

22:28 

15:28 

12.6  ±  1.0 

20:25  i  12 

20:45 

103 

-20  i 

12 

3 

79-04-21 

21:40 

14:51 

12.1  ±  1.5 

19:18  ±  21 

19:21 

136 

-  6  x 

21 

1 

79-11-07 

19:41 

12:35 

9.8  ±  1 .3 

16:08  ±  32 

15:51 

196 

1  t  ill 

32 

5 

79-09-26 

20:04 

13:06 

16.8  ±  2.5 

18:01  ±  19 

17:12 

1  t2 

19  ± 

19 

6 

79-10-08 

23:38 

16:32 

50.0  ±  10.5 

23:10  ±  07 

23:24 

14 

-14  ± 

10 

*'/»•:  Pc  5  wave  onset  time  observed  by  COES  2  given  in  both  universal 
time  ( CT) 

and  magnetic  local  time  (ML'l  ) 

i 7),-.  Estimated  substorm  onset  time  with  uncertainty  i n  units  oi  minutes 
1  /’s :  Actual  substorm  onset  time 
T,;:  Croup  velocity 

All  times  given  to  tens  of  minutes  accuracy. 

Hp  component.  'The  wave  onset  time  at  GOES  3  is  marked  by  an  arrow  in 
Figure  1.  In  the  GOES  2  data,  the  first  peak  of  the  wave  train  that  ran  be 
identified  is  at  21:54  UT  (marked  bv  an  arrow  in  Figure  2).  However,  the 

second  peak  at  21:58  I’T  is  more  pronounced  and  might  be  considered  as  the 

wave  onset.  In  this  case,  the  error  in  determining  the  wave  onset  time  at 
GOES  2  might  be  four  minutes.  Using  21:54  UT  as  the  wave  onset  time  at 
GOES  2.  we  deduced  the  wave  propagation  velocity  to  be  14  km/s.  When  we 
used  21:58  UT  as  the  onset  time  instead,  we  obtained  IS  km/s  for  the  wave 
propagation  velocity.  \Vre  therefore  chose  the  average  propagation  velocity  of 
16  km/s  to  estimate  the  substorm  onset  time,  which  was  found  to  be  20:11 
UT  with  an  uncertainty  of  ±  15  minutes  (case  1  in  Table  1). 

Figure  3  shows  ground  magnetometer  data  of  the  substorm  activity  that, 
might  be  associated  with  the  Pc  5  wave  event  shown  in  Figure  2  (case  l  in 

Table  1).  We  show  only  the  II  component,  of  the  magnetometer  data  from 

five  stations  during  a  14  hour  period  starting  at  17:00  UT.  arranged  from  top 
to  bottom  according  to  geographic  longitude  ranging  from  74°  to  41°  East. 
Near  local  midnight,  indicated  by  an  AM.”  each  station  showed  the  negative 
excursion  characteristic  of  a  substorm.  The  vertical  dashed  line  at  19:55  FT  in 
Figure  3  marks  the  substorm  onset  time  as  determined  by  the  commencement 
of  the  negative  excursion  in  the  II  component  near  midnight.  Note  that  the 
estimated  substorm  onset  time  (20:11  IT)  is  16  minutes  later  than  the  actual 
onset  time  (within  an  uncertainty  of  ±  15  minutes).  Before  this  onset  time, 
ground  magnetometers  show  no  major  substorm  activity.  Therefore.  GOES  2 


GOES  MAGNETOMETER  DATA 


UNI 


Figure  2:  Magnetic  field  measureiT 
dual  Coes  satellites  COES  2  and  3. 


MARCH  25,  1979 


observed  the  Pc  5  event  at  about  15  hours  magnetic  local  tune  approximately 
two  hours  after  the  .substonn  occurred  near  local  midnight. 

As  the  second  example  tease  (i  in  lable  !  i,  one  hour  of  (.'OPS  magne¬ 
tometer  data  starting  at  25:50  H'  October  -S.  1979  are  given  in  Figure  F  [  his 
example  is  difficult  to  study  because  ground  magnetogram--  show  large  sub 
storm  activity.  The  Up  and  He  components  show  the  first  peak  of  the  P<  5 
oscillations  starting  at  25:47  FT  at  GOES  5.  At  the  wave  onset.  GOFS  5  was 
at  11:50  MLT  while  GOES  2  was  at  10:52  MET.  According  to  the  He  compo¬ 
nent,  GOES  2  first  observed  the  lirst  wave  peak  at  25:38  FT  with  the  next 
peak  observed  three  minutes  later.  In  the  same  way  as  case  1.  we  deduced 
that  this  Pc  5  wave  event  had  an  average  propagation  speed  of  50  km/s  with 
an  uncertainty  of  ±  10  km/s.  The  propagation  speed  of  this  event  was  much 
higher  than  that  of  case  1.  Front  the  satellite  location  and  the  group  velocity, 
we  estimated  the  substorm  onset  time  to  be  23:10  FT  with  an  uncertainty  of 
±  7  minutes.  Figure  5  displays  the  ground  data  for  four  stations  in  the  same 
format  as  Figure  5.  The  first  substonn  onset  appeared  to  occur  at  20:40  FT. 
when  the  first  negative  excursion  bay  was  detected  simultaneously  at  three 
stations.  Another  substorm  onset  was  recorded  by  Domhas  station  at  23:24 
FT.  which  is  14  minutes  after  the  estimated  substorm  onset  time.  Therefore, 
the  Pc  5  wave  event  observed  by  the  GOES  satellites  appears  to  be  correlated 
with  the  second  substonn  onset.  We  believe  that  this  event  is  not  correlated 
with  the  first  substonn  onset  because  the  GOES  satellites  had  earlier  detected 
another  Pc  5  event  during  22  FT.  much  closer  to  the  first  substorm  onset. 

Table  1  summarizes  the  results  obtained  by  analyzing  the  six  cases  for 
which  substorm  onset  times  were  reliably  identified,  Comparison  of  the  wave 
onset  time  (given  for  GOES  2  in  Table  1)  with  the  actual  substonn  time 
shows  that  the  wave  events  occurred  generally  2-4  hours  after  the  substonn 
onsets.  However,  one  event  (case  6)  with  a  high  propagation  velocity  (50 
km/s)  occurred  only  14  minutes  after  a  substonn  onset.  Table  1  indicates 
that  a  substorm  onset  occurred  within  20  minutes  of  the  estimated  substonn 
onset  times.  These  agreements  indicate  that  the  Pc  5  waves  and  substonn 
onsets  may  be  correlated. 
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GROUND  MAGNETOMETERS 
H  COMPONENT 


(N.  LAT.,  E.  LONG.) 
MARCH  25-26,  1979 


Figure  3:  Ground  magnetometer  data  in  the  auroral  zone  during  the  stormtime 
Pc  5  event  shown  in  Figure  2.  Above  each  magnetogram  is  the  name  of  the 
station  with  its  geographic  latitude  and  longitude  given  in  parentheses. 
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Figure  4:  Dual  Goes  satellite  observation  of  a  stormtime  Pc  5  wave  event. 
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Figure  5:  Ground  magnetometer  data  in  the  auroral  zone  during  the  stormtime 
Pc  5  event  shown  in  Figure  -4. 
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Ill  Statistical  Survey  of  Propagation  Veloc¬ 
ity  of  Stormtime  Pc  5  waves 

We  have  analyzed  the  GOES  magnetic  field  data  obtained  during  the  years  of 
1979.  1983.  and  1986.  The  GOES  magnetic  field  data  during  these  t  hree  years 
are  selected  because  two  GOES  satellites  were  separated  less  than  two  hours 
local  time.  We  have  scanned  the  GOES  magnetic  field  data  in  the  year  of  1988 
for  Pc  5  magnetic  pulsations  but  did  not  include  the  wave  events  in  the  study 
because  two  GOES  satellites  were  separated  by  more  than  three  hours  local 
time.  When  the  satellite  separation  was  more  than  3  hours  local  time,  two 
GOES  satellites  did  not  detect  Pc  5  wave  events  simultaneously:  often  one 
GOES  satellite  detected  a  clear  Pc  5  event,  while  the  other  GOES  satellite 
merely  detected  a  weak  magnetic  field  depression.  This  is  consistent  with  the 
result  from  the  statistical  study  that  stormtime  Pc  5  wave  events  at  a  single 
observing  point  have  a  duration  of  less  than  2  hours. 

We  scanned  the  magnetic  field  data  and  found  initially  93  events  simul¬ 
taneously  detected  by  two  GOES  satellites.  This  initial  data  set  was  further 
narrowed  down  to  S2  events  by  eliminating  events  unsuitable  for  statistical 
studies.  We  eliminated  events  with  small  wave  amplitudes,  unclear  onsets, 
or  wave  frequency  higher  than  7  mHz.  Note  that  the  wave  period  of  Pc  5 
pulsations  is  generally  defined  in  the  range  of  2  to  10  minutes.  After  elimi¬ 
nating  these  events,  the  data  set  was  reduced  to  71  events.  These  events  were 
then  used  to  deduce  the  statistical  properties.  Appendix  includes  the  lists  of 
stormtime  Pc  5  events  selected  for  the  study. 

A  Survey  Results 

The  histogram  in  Figure  6  shows  that  the  propagation  velocity  of  stormtime  Pc 
o  waves  is  generally  less  than  45  km/s.  The  histogram  of  propagation  velocity- 
distribution  has  a  peak  at  about  15  km/s.  The  wave  period  of  stormtime  Pc 
5  wave  events  distributes  evenly  between  2  and  5  minutes  (Figure  7).  corre¬ 
sponding  to  a  wave  frequency  in  the  range  between  0.001  and  0.007  Hz.  The 
wave  event  tends  to  occur  with  an  average  ambient  magnetic  field  between 
60  and  120  7  (Figure  8)  and  an  inclination  angle  between  15  and  60  degrees 
(Figure  9).  Although  stormtime  Pc  5  events  can  be  detected  from  noon  to 
midnight,  most  frequently  stormtime  Pc  5  waves  are  detected  at  local  time 
between  14  and  20  hours  (Figure  10). 

Figure  11,  which  is  a  histogram  of  wave  event  duration,  indicates  that 
most  events  have  a  duration  less  than  2  hours.  For  each  event,  we  estimated 
the  propagation  velocity  vg  and  the  event  duration  D.  From  vg  and  D .  we 
then  deduced  the  longitudinal  extent  of  wave  region  L  at  synchronous  orbit 
according  to  the  formula  L  -  vgD.  Figure  12  shows  that  a  stormtime  Pc 
5  event  typically  has  a  longitudinal  extent  varying  from  30  to  90  degrees. 
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Figure  6:  Histogram  of  propagation  velocity  of  stormtime  Pc  5  waves. 
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Figure  7:  Histogram  of  wave  period  of  stormtime  Pc  5  waves. 
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Figure  9:  Histogram  of  average  magnetic  field  inclination  angle  during  storm 
time  Pc  5  events. 
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Figure  10:  Distribution  of  stormtime  Pc  5  events  in  local  time. 
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Figure  11:  Histogram  of  stormtime  Pc  5  event  duration. 

However,  a  few  events  have  a  longitudinal  extent  as  large  as  180  degrees. 

According  to  the  statistical  survey,  GOES  satellites  often  observed  storm- 
time  Pc  5  events  in  the  winter  months  from  September  to  April  (Figure  Id). 
The  data  base  contains  no  event  in  June  or  July.  Since  GOES  satellites  were 
stationed  at  the  geographic  equatorial  plane.  GOES  satellites  were  at  high 
geomagnetic  latitudes  during  summer  months.  This  figure  therefore  suggests 
that  stormtime  Pc  5  waves  are  confined  at  low  geomagnetic  latitudes.  This 
conclusion  is  consistent  with  earlier  studies  that  storm  time  Pc  5  waves  have 
an  eigenmode  structure  near  the  equator  (Takahashi  et  al..  19ST).  Theoretical 
studies  presented  in  Section  V  indicate  that  the  wave  mode  is  confined  within 
about  10  degrees  around  the  geomagnetic  equator  near  synchronous  orbit. 

During  stormtime  Pc  5  wave  events,  magnetic  field  magnitudes  are  usually 
comparable  at  two  GOES  satellites  when  they  are  separated  by  two  hours. 
However,  the  inclination  angle  often  varies  drastically  from  the  first  satellite 
(closer  to  midnight)  to  the  second  satellite.  The  histogram  given  in  Figure  14 
indicates  that  the  inclination  angle  at  the  first  satellite  distributed  evenly 
between  Id  and  55  degress,  while  the  inclination  angle  at  the  second  satellite 
was  generally  between  5  and  25  degrees.  This  means  that  the  first  GOES 
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Figure  12:  Distribution  of  longitudinal  extent  for  stormtirr.e  Pc  5  events. 
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Figure  13:  Diurnal  Distribution  of  stormtime  Pc  5  events. 
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Figure  1 4:  Distribution  of  inclination  angle  observed  by  GOES  satellites  during 
stormtime  Pc  5  events.  The  second  GOES  satellite  was  about  30°  west  of  the 
first  GOES  satellite. 

satellite  often  detected  Pc  5  waves  with  a  large  tailward  magnetic  Held,  while 
the  second  GOES  satellite  would  detect  the  event  with  a  more  dipole  magnetic 
field.  This  result  also  implies  that  magnetic  field  configuration  varies  sharply 
Irom  a  strong  tail  field  configuration  to  a  relaxed  dipoie  configuration  within 
a  two  hour  local  time  near  dusk. 

B  Correlation  Study 

Using  the  data  base  and  the  deduced  propagation  velocity,  we  conducted  a  cor¬ 
relation  study  of  propagation  velocity  with  other  observed  parameters.  The 
purpose  of  the  correlation  study  is  to  find  the  dependence  of  propagation  ve¬ 
locity  on  wave  and  ambient  parameters.  The  variation  of  propagation  velocity 
on  parameters  was  then  used  to  determine  the  wave  mode  responsible  for 
stormtime  Pc  5  waves. 

Figure  15  plots  the  propagation  velocity  in  km/s  versus  wave  frequency  in 
Hz.  This  figure  suggests  a  crude  relationship  between  the  propagation  velocity 
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and  wave  frequency:  it  shows  that  the  propagation  velocity  increases  linearly 
with  wave  frequency.  The  correlation  does  not  appear  to  be  strong  because  of 
a  group  of  outer  data  points  with  high  propagation  velocity.  This  may  mean 
that  more  than  one  class  of  wave  modes  are  included  in  the  data  set. 

Furthermore,  Figure  15  shows  that  the  wave  event  scatters  over  a  wide 
range  of  propagation  velocity  for  a  given  wave  frequency.  The  scattering  im¬ 
plies  that  the  propagation  velocity  depends  on  many  parameters.  From  the 
dispersion  equation,  we  expect  that  the  propagation  velocity  depends  on  fre¬ 
quency.  wave  number,  magnetic  field,  plasma  density  and  temperature,  plasma 
pressure,  plasma  beta,  magnetic  field  gradient  and  curvature,  and  other  pa¬ 
rameters. 

To  further  understand  the  dependence  of  propagation  velocity  on  other 
parameters,  we  correlate  propagation  velocity  with  magnetic  field  in  Figure  16. 
which  shows  no  relationship  between  propag.  tion  velocity  and  magnetic  field 
magnitude.  From  this  figure,  we  argue  that  the  propagation  velocity  would 
be  independent  of  magnetic  field  magnitude.  We  therefore  select  events  in  a 
narrow  range  of  magnetic  field  between  70  and  90  7  for  further  studying  the 
coirelation  between  propagation  velocity  and  wave  frequency.  In  this  case,  we 
indeed  obtained  a  much  better  correlation  between  the  propagation  velocity 
and  wave  frequency,  as  shown  in  Figure  17. 

Similarly  when  we  used  the  subset  of  the  data  base  for  7O7  <  B  <  9O7. 
we  found  that  the  propagation  velocity  increases  roughly  with  the  inclina¬ 
tion  angle.  Such  a  correlation  could  not  be  deduced  from  the  complete  data 
set,  probably  for  the  same  reason  that  propagation  velocity  varies  with  many 
parameters. 

Figure  IS  indicates  that  wave  frequency  increases  with  inclination  angle  /. 
which  is  defined  as  tan ~1{D/H).  The  D  and  H  components  are,  respectively, 
the  radial  and  z  component  of  magnetic  fieid  at  the  geographic  equatorial 
plane.  Since  D  is  related  to  the  tail  field,  the  inclination  angle  increases  with 
tail  field.  The  inclination  angle  therefore  measures  the  stretch  of  magnetic 
field  lines  due  to  the  enhancement  of  tail  fields  during  substorms. 

We  have  found  that  wave  frequency  of  stormtime  Pc  5  waves  decreases 
as  the  events  are  detected  closer  toward  noon  (Figure  19).  This  correlation  is 
more  difficult  to  interpret  because  many  parameters  depend  on  local  time.  For 
example,  the  inclination  angle  decreases  from  midnight  toward  neon.  Since 
wave  frequency  is  proportional  to  the  inclination  angle,  wave  frequency  then 
decreases  toward  noon.  However,  the  correlation  between  wave  frequency  and 
local  time  could  be  interpreted  in  several  other  ways.  One  plausible  expla¬ 
nation  is  that  ring  current  ions  injected  during  magnetic  storms  drift  toward 
noon  with  lower  mean  energy  and  a  smaller  density  gradient.  As  a  result,  ring 
current  ions  excite  lower  frequency  waves  toward  noon.  Since  the  plasma  data 
is  not  available  for  study,  we  could  not  make  conclusive  statements  about  this 
possibility. 
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GOES  Pc  5  Wave  Events 


Figure  19:  Wave  frequency  of  stormtime  Pc  5  wave  events  versus  local  time 


IV  Theoretical  Study  of  Low  Frequency  Wave 
Propagation  Velocity 

In  order  to  better  diagnose  substorm  onset  times,  we  have  modeled  the  Pc  5 
wave  propagation  velocity  in  the  equatorial  region.  We  first  solved  the  disper¬ 
sion  equation  of  compressional  waves  derived  from  local  theory.  The  numerical 
results  indicate  that  the  propagation  speed  agrees  with  the  GOES  satellite  ob¬ 
servations.  We  then  compared  the  wave  group  velocities  for  ion  drift  mode, 
drift  mirror  mode,  and  the  low  phase  velocity  drift  mode..  In  addition,  we  in¬ 
vestigated  how  the  group  velocity  varies  with  parallel  and  perpendicular  wave 
lengths. 


A  Dispersion  Equation 

To  model  wave  propagation  velocity,  we  used  a  general  dispersion  equation  for 
an  inhomogeneous  plasma  with  temperature  gradient  in  a  nonuniform  mag¬ 
netic  field  previously  derived  by  Ng  and  Patel  (1983).  This  dispersion  equation 
has  previously  been  used  to  study  drift  wave  instabilities  in  the  magnetosphere. 

The  plasma  is  assumed  to  have  a  cold  plasma  density  nc  and  a  hot  plasma 
density  n/,,  a  density  gradient  scale  length  rn,  a  magnetic  field  gradient  scale 
length  r\i,,  plasma  j3,  and  temperature  anisotropy  a.  We  assume  the  ratio 
between  the  density  gradient  scale  length  and  the  perpendicular  temperature 
gradient  scale  length  to  be  77.  The  dispersion  equation  is  then 

D\\D27  —  ~^-DuD2i  =  0  (2) 


where 


D 


11 


4. «. 


% 

+—(T0-rfb{TQ-ri))-~ 

l o  uj 


(i_^(i_r?)j(r0-6(r0-rl)) 


ijj. 


D 12  — 


u>2 


+  -^1; 


'-r,  +  a—  )  ((2  -  b)V0  -  6(3  -  26) (T0  -  t\)) 

"  1  \  OJ  /  \  U)  uJ 


nh<U 


(To  +  (i  -  26)(r0  -  r,))  -  —l2j\ 

UJ  J 


29 


k2 


^22  —  ^2  y7  ~  T.  $1.1  MTo  —  T\)  +  L$ 


The  parameters  are  defined  as 
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7\  and  7jj  are  respectively  the  perpendicular  and  parallel  temperatures,  Z  is 
the  usual  plasma  dispersion  function  with  argument  (u>  ~ ui^t) / {k^\ j^),  J0  and 
Ji  are  the  zeroth  and  first-order  Bessel  functions  with  arguments  k^ajt1^2;  /0 
and  I\  are  modified  Bessel  functions  with  arguments  bj.  The  drift  frequencies 
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B  Numerical  Solutions 

For  a  given  parallel  wave  length,  we  first  solved  numerically  Equation  2  for 
wave  frequency  as  a  function  of  perpendicular  wave  number  kx.  We  then 
calculated  the  perpendicular  group  velocity  Vg  from  the  relationship 

Vg  =  duj/dkx  (3) 

The  propagation  velocity  is  examined  for  three  wave  modes:  ion  drift,  drift 
mirror  and  the  low  phase  velocity  drift  modes.  The  parameters  used  in  the 
calculation  are:  density  ratios  nhi/nct  and  nhe/nce  =  0.01,  t]i  and  7je  =  —3.0, 
rn  —  50  and  temperature  anisotropy  A  =  0,  @\\  =  1,  T,  =  10  keV,  and 
bt  =  l /2(kxpl)2  =  0.1.  The  solutions  of  the  dispersion  equation  for  various 
parameters  are  shown  in  Figures  20-26. 

B.l  Ion  Drift  Mode 

Figure  20  shows  the  wave  frequency  (solid  line)  and  growth  rate  (dashed  line) 
as  a  function  of  perpendicular  wave  vector  times  the  ion  gyroradius  kxp,  for 
the  ion  drift  mode.  For  the  ion  drift  mode,  the  frequency  is  positive  for  positive 
wave  mode  number  and  thus  the  wave  phase  velocity  is  in  the  direction  of  the 
ion  diamagnetic  drift.  For  this  mode,  the  parallel  wave  number  is  very  small 
and  thus  the  wave  mode  has  a  long  parallel  wave  length.  This  figure  indicates 
that  the  instability  occurs  at  very  small  parallel  wave  number  (A'up,  <  0.006) 
and  the  growth  rate  decreases  with  The  wave  frequency  of  growing  waves 
u>  normalized  by  A.'||U,  decreases  from  3  to  1  as  k^pi  increases  from  0.002  to 
0.006. 
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Figure  21:  Perpendicular  group  velocity  as  a  function  of  the  parallel  wave 
vector  for  the  ion  drift  mode. 

However,  according  to  the  GOES  observations,  the  wave  is  confined  to 
10°  latitude.  The  parallel  wave  number  A’np,  should  be  greater  than  0.04  for  a 
typical  ion  temperature  of  10  keV.  Therefore,  ion  drift  waves  have  difficulties 
in  explaining  the  parallel  wavelengths  of  the  Pc  5  waves  observed  by  GOES 
satellites. 

For  the  ion  drift  mode,  the  group  velocity  ir'cr-=>3ce«  *he  parallel  wave 
number  fc||  (Fig  21).  For  A||pj  =  0.002.  when  the  ion  drift  mode  has  a  large 
growth  rate,  the  perpendicular  group  velocity  is  less  than  10  km/s.  Since  most 
propagation  velocity  for  stormtime  Pc  5  waves  are  greater  than  10  km/s,  the 
GOES  satellite  observations  suggest  that  ion  drift  mode  could  not  account  for 
stormtime  Pc  5  waves  with  large  propagation  velocity. 

In  Figure  22,  we  plot  the  frequency  and  growth  rate  as  a  function  of  kLpi 
for  the  ion  drift  mode.  The  parameters  are  the  same  as  Figure  20  except 
that  k\\pi  =  0.002  is  chosen.  This  figure  indicates  that  the  ion  drift  waves  are 
unstable  for  k±p,  <  1.  In  Figure  23.  we  next  present  the  perpendicular  group 
velocity  as  a  function  of  k±p,  for  the  ion  drift  mode.  The  perpendicular  group 
velocity  increases  from  10  km/s  to  about  20  km/s  as  kj_p,  increases  from  0.1 
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Figure  22:  Wave  frequency  and  growth  rate  as  a  function  of  the  perpendicular 
wave  vector  for  the  ion  drift  mode. 
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Figure  23:  Perpendicular  group  velocity  as  a  function  of  the  perpendicular 
wave  vector  for  the  ion  drift  mode. 
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Figure  24:  Wave  frequency  and  growth  rate  as  a  function  of  the  parallel  wave 
vector  for  the  drift  mirror  mode. 

to  1.  Since  Pc  5  waves  generally  have  a  wave  propagation  velocity  less  than  50 
km/s.  it  is  likely  that  the  perpendicular  wavelength  needs  to  be  greater  than 
the  ion  gyroradius  such  that  k^p,  <1.  It  is  generally  expected  that  kp,  is  less 
than  1.  The  statistical  survey  suggests  that  the  perpepndicular  propagation 
velocity  varies  from  5  km/s  to  40  km/s  and  for  about  half  of  the  event'  the 
propagation  velocity  is  less  than  20  km/s. 

B.2  Drift  Mirror  Mode 

We  next  investigate  propagation  velocity  for  the  drift  mirror  mode.  Figure  21 
plots  the  frequency  and  growth  rate  as  a  function  of  k^p,  for  the  drift  mirror 
mode.  For  the  drift  mirror  mode,  which  is  unstable  for  a  plasma  with  a 
temperature  anisotropy,  we  chose  the  following  parameters:  density  ratios 
and  iuir/  n.,,  =  0.01.  density  gradient  scale  length  r„  =  50,  temperature 
anisotropy  A  =  1.5,  Jj|  =  1.  T,  =  10  keV.  and  b,  =  1/2 {kLp,)2  —  0.1.  To 
separate  the  effects  of  drift  mirror  mode  from  those  of  ion  drift  mode,  we 
examie  the  case  in  the  absence  of  temperature  gradients  ( r;,  and  ?/,  =  0). 
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Figure  25:  Perpendicular  group  velocity  as  a  function  of  the  parallel  wave 
vector  for  the  drift  mirror  mode. 

For  a  given  kLpt,  the  drift  mirror  waves  are  unstable  for  all  Arjjp,,  and  the 
wave  frequency  is  much  less  than  k^vx  {  ^ / k\\ vx  <<  1)  (Figure  211.  The  growth 
rate  y/k^Vi  also  decreases  with  increasing  For  the  same  parameters,  the 

perpendicular  group  velocity  as  a  function  of  k^p,  for  the  drift  mirror  mode  is 
shown  in  Figure  25.  As  k||/?i  increases  from  0.05  to  0.3.  the  perpendicular  group 
velocity  decreases  from  60  km/s  to  about  10  km/s.  In  the  range  of  k^px  <  0.5. 
the  perpendicur  group  vlocitv  is  between  5  to  60  km/s.  The  calculted  range 
of  perpendicular  group  velocity  therefore  agrees  with  the  observations. 

We  next  examine  the  dependence  of  drift  mirror  mode  on  wave  perpendic¬ 
ular  wavelength.  Figure  26  plots  the  frequency  and  growth  rate  as  a  function 
of  kipi  for  the  drift  mirror  mode.  In  this  figure,  we  used  k\\px  —  0.2  and  other 
parameters  the  same  as  Figure  24.  Figure  27  plots  the  perpendicular  group 
velocity  as  a  function  of  kLpx  for  the  drift  mirror  mode.  The  perpendicular 
group  velocity  increases  from  5  km/s  to  about  30  km/s  as  k^_p,  increases  from 
0.1  to  0.6.  In  the  range  of  kLpx  <  l,  the  perpendicur  group  vlocitv  li,i  is 
between  5  to  30  km/s.  The  perpendicular  group  velocity  i^  maximum  when 
k±p,  —  0.6.  For  kLpi  >  0.6.  perpendicular  group  velocity  V3±  decreases  with 
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Figure  28:  Wave  frequency  and  growth  rate  as  a  function  of  the  ratio  between 
perpendicular  and  parallel  wave  vectors  for  the  drift  mirror  mode. 

increasing  k± p,. 

The  drift  mirror  mode,  which  is  excited  by  temperature  anisotropy,  has 
the  largest  growth  rate  at  small  A.'|j / (Figure  2S).  Figure  29  shows  that  the 
group  velocity  of  the  drift  mirror  mode  has  a  wide  range,  varying  from  about 
50  km/s  at  small  k^/k^  to  zero  when  k\\/ki_  is  greater  than  1.  Therefore  the 
drift  mirror  mode  agrees  well  with  the  observed  range  of  propagation  velocity 
for  stormtime  Pc  5  waves,  which  varies  from  from  5  km/s  to  40  km/s. 


B.3  Low  Phase  Velocity  Drfit  Mode 

The  temperature  gradient  can  also  excite  a  low  frequency  wave  mode  with  low 
phase  velocity.  This  low  phase  velocity  mode  differs  from  the  ion  drift  mode 
in  having  a  transverse  magnetic  field  component.  For  this  mode,  the  value  of 
transverse  to  compressional  magnetic  field  amplitudes  reaches  as  high  as  0.6. 
whereas  the  ion  drift  mode  is  mainly  compressional. 

Figure  30  indicates  that  the  instability  of  low  phase  velocity  drift  wave 
mode  occurs  for  all  values  of  k±p,.  When  k±p,  <  1,  the  wave  frequency  is 
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Figure  29:  Perpendicular  group  velocity  as  a  function  of  the  ratio  between 
perpendicular  and  parallel  wave  vectors  for  the  drift  mirror  mode. 
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Figure  30:  Wave  frequency  and  growth  rate  as  a  function  of  the  perpendicular 
wave  vector  for  the  low  phase  velocity  drift  mode. 

positive.  The  wave  frequency  becomes  negative  as  kj_px  increases  above  1. 
Figure  31  shows  that  the  group  velocity  for  this  low  phase  velocity  mode  is 
generally  less  than  10  km/s.  Furthermore,  the  group  velocity  is  positive  only 
when  k±pt  <  0.7  (positive  velocity  corresponds  to  westward  propagation  in  our 
geometry).  Since  storm  time  Pc  5  waves  generally  propagate  westward  with 
a  velocity  greater  than  10  km/s.  this  wave  mode  cannot  satisfactorily  explain 
the  GOES  satellite  observations. 
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V  Eigenmode  Analysis  of  Low  Frequency  Waves 
in  a  Dipole  Magnetic  Field 


A  Eigenmode  Equation 

When  the  parallel  wavelength  is  comparable  to  the  curvature  scale  length,  the 
local  dispersion  equation  is  no  longei  a  good  approximation.  We  therefore 
proceeded  to  solve  an  eigenmode  equation  of  compressional  waves  in  a  dipole 
magnetic  field.  We  used  a  gyrokinetic  formalism  to  derive  the  eigenmode 
equation  (Cheng  and  Lin,  1987) 


B  •  V 
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2F2 
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x[cos(p6)(co$(pQ)B\\)  -f  sin(p0)(sin(p0)B||)] 
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=  0 


(4) 


Briefly  we  list  the  definitions  in  the  equation: 

(1)  By  is  the  parallel  component  of  the  complex  wave  field  along  ambient  field 
B  normalized  by  the  magnitude  of  the  equatorial  ambient  field. 

(2)  u>  is  the  complex  eigenfrequency. 

(3)  F  =  (u;  -  (w))/u.v 

(4)  (. . .)  denotes  an  averaging  quantity  along  a  trapped  ion  trajectory  for  a 
one  bounce  cycle. 

(5)  0  is  a  function  of  polar  angle  and  is  expressed  as  an  integral  over  a  function 
of  the  ambient  field. 

The  first  two  terms  of  Equation  4  give  the  conventional  expression  of  com¬ 
pressional  waves  from  magnetohydrodynamic  (MHD)  theory  without  kinetic 
effects.  The  integral  contains  the  kinetic  effects  contributed  by  the  trapped 
ions.  This  eigenmode  equation,  which  can  be  solved  for  the  eigenfunctions 
of  B||  and  the  frequency  u;,  means  physically  that  Pc  5  waves  can  be  ap¬ 
proximated  as  standing  waves  in  the  magnetosphere  along  ambient  field  lines. 
Overall,  Equation  4  is  an  integro- differential  equation  that  requires  a  five¬ 
dimensional  integration  (3  velocity  integrations,  1  bounce  average  integration 
and  the  integration  for  0).  Furthermore,  unlike  the  local  dispersion  equation, 
the  frequency  appears  in  a  highly  non-linear  way  through  the  function  F.  All 
this  makes  solving  the  equation  computationally  intensive.  If  we  neglect  the 
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last  term  with  the  integral  (—  -iir  fT  ■  ■  ■)  which  takes  into  account  kinetic  effects, 
the  problem  is  simplified  to  the  MHD  case: 
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Our  model  uses  a  dipole  magnetic  field 


B  =  —(2  sin  A  f  —  cos  A  X  ) 


(5) 


(6) 


On  a  given  held  hne,  r  —  L  cos2  A,  and  we  introduce  the  variable  x  =  sin  A,  so 
that  our  differential  equation  can  be  written  as 

<l(x)-J—2+  p{x)~r  ~  Hz)  -  u;2s(x)]  =  G,  (7) 


Here  we  have  defined  b(x)  —  £?j|(A).  Other  coefficients  in  Equation  7  are 
defined  as: 
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We  also  use  the  following  definitions: 

<7  =  1+  (&x  +  0i\)f  2 

h{x)  =  BojB{x) 


kj_  perpendicular  wave  number 

L  geosynchronous  orbit  radius  («  6.6 Re) 

$U  parallel  fluid-to-magnetic  pressure  ratio 

0x  perpendicular  fluid-to-magnetic  pressure  ratio 

Ti|  parallel  temperature 

TL  perpendicular  temperature 

uy,  ion  cyclotron  frequency 


Alfven  speed. 


v;4 

The  boundary  conditions  are  chosen  at  x  =  0  and  at  x  =  x\t  such  that 

b(  0)  =  0 

b(xM)  =  0 


In  the  numerical  solution  we  use  xM  —  sin  20°. 

To  solve  the  equation  we  use  the  finite  difference  method.  We  solve  for  6 
values  at  N  evenly-spaced  points  between  0  and  xM,  so  that  Ax  =  xA,/(iV  -f 1 ). 
Specifically,  we  approximate  the  derivatives  at  each  point  by  central  differ¬ 
ences: 
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(Ax)3 


where  the  subscript  i  means  the  value  at  x,  =  i  Ax.  We  take  the  finite  difference 
of  Equation  7  at  each  point  x;  to  get  a  system  of  N  equations  in  the  (N  +  1 ) 
unknowns  b^,oj: 


,  P*  )  L  . 
(Ax)2  +  2Ax  '  1+1  + 


2  q, 


(Ax)2 


~r  —  r,  -f-  J2sl  )  6,  + 


(Ax)2 


or 


b,-i  -  0 


for  i  =  1 _ ,  N  and  where  b0  =  0  and  V+i  =  0  from  our  boundary  conditions. 

In  matrix  notation  we  have 


A(w)  •  b  =  0  (8) 

in  which  A(u>)  is  an  N  x  N  tridiagonal  matrix.  If  b  has  nontrivial  solutions, 
we  obtain 


det[A(w)J  =  0.  (9) 

We  numerically  solve  this  equation,  using  Newton  method  for  solving  the  roots 
to  obtain  eigenfrequencies  u>.  For  a  particular  solution  w,  the  system  of  equa¬ 
tions  is  linear  and  homogeneous,  so  we  expect  to  determine  the  eigenfunction 
up  to  a  constant  factor.  We  normalize  our  solutions  by  taking  6j  =  l.  Using 


this  approach,  we  find  a  range  of  eigenfrequencies  similar  to  those  obtained  by 
analytical  solution  of  a  local  approximation  to  the  MHD  equation. 

We  next  describe  the  evaluation  of  the  kinetic  term 
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We  first  rewrite  the  integral,  using  the  relation 
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where  A  =  h(xT),  k  at  the  trapped-  particle  turning  point.  We  next  substitute 
the  variables  Fh  and  7j|  in  terms  of  v  =  v/v,  A  =  (7x/7j|)o,  and  /?ii  according 
to  the  following  definitions: 
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where  the  subscript  0  means  to  evaluate  at  zero  latitude. 

This  substitution  allows  us  to  express  the  kinetic  term  in  terms  of  known 
quantities.  In  the  case  of  odd  modes,  when  (B\\)  =  0,  we  express  the  kinetic 
term  as 
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L'o  derive  the  computational  formulas  in  terms  of  i\  we  use  the  following 
definitions: 
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and  the  following  relationships: 
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Using  these  definitions,  we  derive  the  variables  in  terms  of  x: 
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where  the  over  bar  means  that  the  frequency  has  been  scaled  by  and  Lh 
is  the  pressure  gradient  scale  length.  To  compute  the  bounce  averages  of  a 
variable  /,  we  use  the  definition 


rt::iTdi/v  „ 


Since  B\\  is  an  odd  function  of  x ,  cos(p8)B\\  is  also,  and  we  have 


W  =  0 


(cos(p0)Bj|)  =  0. 

On  the  other  hand,  sin(p8)Bji  and  uj  are  even  functions  of  x.  and  we  evaluate 
them  using 


</)  -  ZSMptl. 

In  evaluating  the  integral  /(d//u||)  we  encounter  a  singularity  at  x  =  xT. 
Near  this  point  we  use  an  analytic  function  to  approximate  the  integrand  and 
integrate  the  approximation  in  closed  form. 

When  computing  (s'm{p6)B\\)  the  period  of  sin(p(?)  may  be  on  the  or¬ 
der  of  the  finite-difference  spacing.  To  accurately  evaluate  the  integral,  we 
compute  the  integral  over  each  subinterval  from  the  known  values  of  B\\  and 
(!/u|l )(dt/dx)  at  each  of  the  evenly-spaced  points.  We  assume  that  these 
functions  are  linear  over  each  subinterval,  and  we  integrate,  in  closed  form. 
s\n(p0) B\\{1  f  v\\){dl / dx).  Then  we  simply  sum  these  integrals  over  all  subinter- 
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vals. 

B  Numerical  Solutions 

The  procedure  of  solving  the  eigenmode  equation  is  organized  as  follows:  First, 
we  solve  the  MHD  case,  as  outlined  above.  The  solution  is  then  used  as  an 
initial  guess  in  our  iterative  scheme.  We  proceed  to  solve  for  a  new  eigenfre- 
quency,  using  the  kinetic  term  computed  from  the  MHD  eigenfrequencv  and 
eigenfunction.  Then  we  calculate  a  new  eigenfunction,  corresponding  to  the 
new  eigenfrequencv.  This  procedure  is  iterated  until  the  estimated  error  is 
within  acceptable  bounds. 

Soving  the  eiegnmode  equation,  we  obtained  the  wave  frequency  and 
growth  rate  of  the  unstable  eigenmode.  Figure  32  plots  the  real  and  imag¬ 
inary  frequencies  as  a  function  of  anisotropy.  In  addition,  the  growth  rate 
solving  from  the  MHD  equation  is  also  plotted  in  Figure  32  for  comparison. 
In  the  MHD  limit,  the  frequency  of  drift  mirror  waves  is  zero.  Including  the  ki¬ 
netic  effects,  the  drift  mirror  mode  has  a  small  frequency,  about  0.01  of  the  ion 
cyclotron  frequency  (w/wy,  —  0.01).  More  importantly,  the  eigenmode  equatin 
indicates  a  smaller  growth  rate,  in  better  agreement  with  the  observations. 

Figure  33  shows  the  eigenmode  structure  of  magnetic  field  for  the  drift 
mirror  mode  versus  0  for  the  parameters  :  J||  =  0.5. A  =  2.2.  Lh/p  —  —50. 
ki_L  =  64.  rih/nc  =  0.1, L/p  —  425.  7 ]\/Tc  =  1000.  and  Lh/Lh  =  -0.1.  This 
figure  indicates  that  the  wave  amplitude  reaches  maximum  at  0  ~  V'  and 
decreases  to  zero  when  0  >  10°.  The  wave  amplitude  of  drift  mirror  waves  is 
therefore  localized  near  the  equator.  Since  the  wave  amplitude  vanishes  near 
the  equator,  the  wave  has  an  odd  mode  structure  with  respect  to  the  equator. 
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Figure  32:  Wave  frequency  and  growth  rate  as  a  function  of  anisotropy  for  t  he 
drift  mirror  mode 
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Eigenmode  Structure 


Latitude  (degrees) 


Figure  33:  Eigenmocle  structure  of  magnetic  field  oscillations  for  the  drift 
mirror  mode.  The  y  axis  is  the  magnetic  field  wave  amplitude  normalized  bv 
magnetic  field  at  the  equator,  and  the  x  coordinate  is  the  geomagnetic  latitude. 


VI  Discussion 


This  project  has  demonstrated  that  stormtime  Pc  5  waves  observed  by  geosta¬ 
tionary  satellites  in  the  afternoon  sector  are  correlated  with  substorm  onsets 
occurring  near  local  midnight.  The  disturbed  region  of  a  stormtime  Pc  5  event 
has  a  longitudinal  extent  varying  between  30  and  90  degrees.  The  study  shows 
that  stormtime  Pc  5  waves  have  a  wave  amplitude  confined  with  about  10° 
from  the  magnetic  equator.  The  propagation  velocity  of  stormtime  Pc  5  waves 
is  typically  about  15  km/s  and  can  be  as  high  as  45  km/s.  The  numerical  so¬ 
lutions  suggest  that  the  propagation  velocity  of  stormtime  Pc  5  waves  agrees 
better  with  the  perpendicular  group  velocity  of  drift  mirror  mode. 

The  correlation  between  the  occurrence  of  stormtime  Pc  5  wave  events 
and  substorm  onsets  is  not  obvious  because  stormtime  Pc  5  waves  generally 
occur  2-4  hours  after  substorm  onsets.  Using  the  wave  propagation  velocity, 
we  estimated  the  time  for  the  wave  to  propagate  from  local  midnight  to  the 
satellites.  For  the  events  studied,  we  were  able  to  identify  a  substorm  onset 
from  ground  magnetograms  near  local  midnight  within  20  minutes  of  the  esti¬ 
mated  times.  Therefore,  the  results  suggest  that  the  occurrence  of  stormtime 
Pc  5  waves  is  correlated  with  substorm  onsets.  The  correlation  implies  that 
stormtime  Pc  5  waves  are  propagating  westward  with  the  substorm  injected 
plasma.  Since  the  propagation  velocity  is  related  to  the  energy  of  injected 
plasma,  the  larger  propagation  velocity  would  indicate  that  GOES  satellites 
will  encounter  more  energetic  plasma  as  they  move  toward  local  midnight. 

The  correlation  results  depend  on  the  accuracy  in  determining  the  wave 
propagation  velocity.  In  order  for  GOES  satllites  to  detect  accurately  the 
wave  propagation  velocity,  the  separation  between  the  two  GOES  satellites 
needs  to  be  less  than  2  hours.  We  believe  that  the  error  in  deducing  the  wave 
propagation  speed  is  mainly  caused  by  uncertainties  in  identifying  the  time  of 
the  first  oscillation  peak,  which  could  at  most  have  an  error  of  one  wave  period. 
As  Table  1  shows,  the  error  in  propagation  speed  is  generally  a  few  km/s  when 
the  propagation  speed  is  greater  than  10  km/s.  This  produces  uncertainties 
in  estimating  the  substorm  onset  times  by  about  ±  20  minutes  (Table  1). 
Since  the  actual  substorm  onsets  generally  occurred  within  20  minutes  of  the 
estimated  substorm  onset  time,  the  uncertainties  in  the  propagation  speed  do 
not  affect  the  qualitative  conclusion  that  stormtime  Pc  5  waves  are  related  to 
substorm  onsets. 

In  addition  to  uncertainties  in  propagation  speed,  other  factors  may  af¬ 
fect  the  estimate  of  substorm  onset  times.  For  example,  we  have  assumed  a 
constant  wave  group  velocity  in  the  estimate,  but  the  wave  group  velocity  may 
vary  with  longitude  and  time,  depending  on  the  azimuthal  variation  of  plasma 
parameters.  Because  the  GOES  2  and  3  satellites  at  the  geosynchronous  orbit 
are  separated  only  by  two  hours  local  time,  the  effects  of  the  radial  motion  on 
estimating  wave  propagation  speed  should  be  small.  Furthermore,  since  most 
stormtime  Pc  5  waves  occurred  in  the  afternoon  sector  before  dusk,  the  radial 
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motion  of  the  injected  plasma  boundary  may  be  unimportant  for  these  events. 

Another  interesting  observation  is  that  the  GOES  2  satellite  often  detected 
a  stormtime  Pc  5  event  not  observed  by  the  GOES  5  satellite.  This  is  peculiar 
because  stormtime  Pc  5  waves  presumably  propagate  westward  from  GOES  5 
to  GOES  2.  When  GOES  5  did  not  observe  the  Pc  5  events,  the  magnetic 
field  at  GOES  5  had  a  large  tail  field.  The  explanation  might  be  that  the 
GOES  5  satellite  was  located  at  a  higher  magnetic  latitude  than  GOES  2  and 
the  waves  were  localized  near  the  magnetic  equator.  This  could  also  explain 
the  observation  that  stormtime  Pc  5  waves  are  not  detected  at  large  inclination 
angle  (Baefield  and  Lin.  1983:  Lin  and  Cheng.  1984J. 

The  size  of  the  region  of  stormtime  Pc  5  waves  has  not  been  known  be¬ 
fore.  The  information  about  the  size  of  disturbance  region  generally  cannot 
be  obtained  from  single  satellite  observations.  By  using  dual  satellite  obser¬ 
vations  about  the  propagation  velocity  and  event  duration,  we  were  able  to 
estimate  that  the  disturbed  region  of  a  stormtime  Pc  5  event  has  a  longitudi¬ 
nal  extent  varying  between  30  and  90  degrees.  The  wave  disturbance  region 
in  the  evening  hours  is  on  the  average  smaller  the  disturbance  region  in  the 
afternoon  hours.  The  plasma  clouds  injected  during  substorms  are  responsible 
for  exciting  stormtime  Pc  5  waves.  As  the  substorm  injected  plasma  cloud 
drifts  westward,  the  front  moves  faster  than  the  trailing  edge.  As  a  result,  the 
plasma  cloud  expands. 

The  numerical  solutions  of  the  dispersion  equation  suggest  that  the  prop¬ 
agation  velocity  of  stormtime  Pc  5  waves  can  be  explained  by  the  drift  mirror 
mode  (Section  IV).  The  range  of  propagation  velocity  can  be  explained  when 
the  perpendicular  wave  vector  times  the  ion  gyroradius  is  less  than  1  (k±p,  <  1) 
and  the  parallel  wave  vector  times  ion  gyroradius  is  in  the  range  between  0.05 
and  0.3  (0.05  <  k^p,  <  0.3)  (see  Figure  25  and  27). 

The  correlation  study  indicates  that  the  propagation  velocity  increases 
with  wave  frequency.  Based  on  the  numerical  results,  we  argue  that  this 
correlation  can  be  explained  if  parallel  wavelength  is  the  dominant  parameter 
affecting  the  wave  frequency  of  stormtime  Pc  5  waves.  Figures  24  and  25 
indicate  that  both  wave  frequency  and  propagation  velocity  decrease  with 
fc|]/V  Therefore,  propagation  velocity  would  increase  with  wave  frequency  as 
parallel  wavelength  increases.  The  correlation  between  propagation  velocity 
and  wave  frequency  cannot  be  explained  by  varying  perpendicular  wavelength 
according  to  Figures  26  and  27.  These  two  figures  show  that  the  perpendicular 
group  velocity  increases,  whereas  wave  frequency  decreases  as  k±pt  increases. 

The  correlation  would  also  require  that  other  parameters  including  per¬ 
pendicular  wavelength  and  plasma  parameters  do  not  vary  much.  This  imply 
that  the  wave  frequency  of  stormtime  Pc  5  wave  events  is  mainly  determined 
by  parallel  wavelength.  The  eigenmode  analysis  presented  in  Section  V  has 
demonstrated  that  parallel  wavelength  of  low  frequency  waves  in  a  dipole  filed 
is  not  an  independent  parameter.  Parallel  wavelength  is  an  independent  pa¬ 
rameter  only  when  the  ambient  magnetic  field  is  straight.  Parallel  wavelength 
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is  determined  by  the  magnetic  field  configuration  or  magnetic  field  curvature. 

The  correlation  study  between  wave  frequency  and  magnetic  field  proper¬ 
ties  supports  the  eigenmode  analysis.  Since  wave  frequency  of  stormtime  Pc  5 
waves  is  also  correlated  with  the  inclination  angle  (Figures  17  and  IS),  parallel 
wavelength  is  in  turn  determined  by  the  inclination  angle  or  the  magnetic  field 
topology.  A  magnetic  field  configuration  with  a  large  tail  field  during  magnetic 
storms  will  have  a  large  inclination  angle  or  a  large  magnetic  field  curvature. 
For  such  a  magnetic  field  configuration,  stormtime  Pc  5  waves  would  have  a 
higher  wave  frequency. 

In  summary,  the  propagation  velocity  of  stormtime  Pc  o  waves  is  consis¬ 
tent  with  the  perpendicular  group  velocity  of  drift  mirror  mode.  The  propaga¬ 
tion  velocity  increases  with  wave  frequency  and  the  magnetic  field  inclination 
angle.  The  wave  frequency  of  stormtime  Pc  5  waves  appears  to  be  mainly 
determined  by  wave  parallel  wavelength,  which  is  in  turn  determined  by  the 
inclination  angle  or  the  magnetic  field  topology. 

Finally  the  present  study  suggests  that  remote  diagnosis  of  substorm  ac¬ 
tivities  from  dayside  snchronous  satellites  is  feasible.  The  onset  of  stormtime 
Pc  5  waves  observed  in  the  dayside  can  be  used  to  help  specify  the  night- 
side  ionospheric  condition,  and  the  disturbed  environment  at  the  synchronous 
orbit.  In  principle,  it  is  feasible  to  use  real  time  observations  of  magnetic 
pulsations  to  diagnose  remotely  the  plasma  environment  that  a  synchronous 
satellite  is  about  to  go  through  during  substorms.  For  example,  a  fast  moving 
stormtime  Pc  5  event  indicates  that  the  magnetic  field  configuration  ahead  in 
the  evening  sector  is  stretched  tailward.  Therefore,  more  substorm  injections 
are  likely.  On  the  other  hand,  a  slow  moving  stormtime  Pc  5  event  might 
indicate  a  more  relaxed  magnetic  field  configuration.  Thus  substorm  activity 
is  subsided.  The  observations  of  low  frequency  magnetic  pulsations  might  be 
useful  for  planning  synchronous  satellite  operations 
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Appendix 

List  of  Stormtime  Pc  5  Events 


List  of  Stormtime  Pc  5  Wave  Events  in  1979  -  1980 


No. 

Date 

Start  Time 
(HH:MM) 

Time  Delay 
(min) 

Wave 

period 

(min) 

Velocity 

(km/s) 

1 

3/25/79 

21:38 

4.5 

6.67 

2 

3/28/79 

20:18 

9.25 

8.87 

3 

3/28/79 

22:28 

26 

3.0 

11.17 

4 

3/28/79 

23:08 

34 

2.75 

7.82 

5 

3/29/79 

18:00 

17 

6.8 

18.72 

6 

4/01/79 

17:18 

30 

7.4 

9.27 

7 

4/02/79 

16:48 

43 

9.5 

5.53 

8 

4/21/79 

21:37 

31 

3.8 

8.87 

9 

9/20/79 

19:24 

10 

6.1 

34.61 

10 

9/21/79 

0:36 

2.4 

50.59 

11 

9/26/79 

1:24 

2.8 

21.80 

12 

9/26/79 

20:04 

4.8 

18.87 

13 

9/29/79 

0:10 

3.5 

10.24 

14 

10/8/79 

23:38 

3.0 

37.98 

15 

10/12/79 

23:06 

24 

4.4 

12.21 

16 

11/13/79 

20:03 

19 

5.4 

15.32 

17 

11/15/79 

18:35 

1 1 

9.2 

28.82 

18 

11/24/79 

21:42 

14 

4.0 

22.07 

19 

1/1/80 

19:03 

42 

5.6 

4.94 

20 

1/1/80 

21 :42 

25 

4.4 

10.41 

21 

1/13/80 

18:54 

25 

6.0 

10.31 

22 

1/13/80 

22:03 

21 

3.5 

12.86 

23 

1/14/80 

1:18 

20 

4.4 

13.60 

24 

2/16/80 

19:34 

37 

6.0 

5.66 

*  Event  date  and  start  time  were  based  on  GOES  2  observations. 

*  time  delay  is  measured  from  the  wave  onset  time  at  GOES  2  to  the 
wave  onset  time  at  GOES  3. 


List  of  Stormtime  Pc  5  Wave  Events  in  1983 


Date 

Start  Time 

Time  Delay 
(min) 

1/14/83 

23:59 

25 

1/15/83 

01:51.5 

21.5 

1/17/83 

23:21 

12 

1/24/8J 

22:46.5 

13 

2/11/83 

20:25 

13 

2/20/83 

19:44 

57.5 

3/5/83 

22:13 

21 

3/12/83 

18:26 

41 

3/18/83 

22:22.5 

19.5 

3/25/83 

20:03 

39 

3/28/83 

17:19.5 

10 

3/28/83 

19:15 

9 

4/14/83 

20:12 

10.5 

4/14/83 

23:44.5 

15.9 

4/24/83 

23:07.5 

10.5 

4/30/83 

21:35 

17 

9/19/83 

1:03 

22 

10/14/83 

1:24 

26 

10/18/83 

20:06 

24 

10/24/83 

0:18 

20 

11/12/83 

20:28 

34 

11/15/83 

12:29 

*28 

11/25/83 

22:33 

41 

12/5/83 

22:19 

48 

12/31/83 

23:00 

22 

3/1/83 

0:13 

28 

3/2/83 

14:36 

57 

3/12/83 

20:07 

62 

3/28/83 

21:07 

48 

13.52 

16.23 

31.52 
28.86 
28.38 

3.98 

16.10 

6.71 

17.70 

7.31 

37.58 

42.10 

35.64 

22.49 

35.52 
20.61 
30.40 

25.53 
28.02 
34.36 
19.12 

-30.02 

14.94 

12.23 
31.02 
11.30 

3.98 

3.39 

5.38 


Event  date  and  start  time  were  based  on  GOES  5  observations, 
time  delay  is  measured  from  the  wave  onset  time  at  GOES  5  to  the 
wave  onset  time  at  GOES  2. 


List  of  Stormtime  Pc  5  Wave  Events  in  1986 


No. 

Date 

■ 

Time  Delay 
(min) 

Wave 

period 

(min) 

1 

1/6/86 

20:59 

47 

4.0 

5.67 

2 

1/9/86 

22:00 

25 

3.0 

13.47 

3 

1/21/86 

22:11 

56 

3.4 

4.37 

4 

1/25/86 

3:24 

15 

8.5 

24.85 

5 

2/11/86 

4:10 

32 

2.7 

10.04 

6 

2/21/86 

23:06 

13 

3.0 

29.14 

7 

3/7/86 

20:31 

15 

3.2 

24.35 

8 

3/25/86 

22:16 

15 

3.6 

23.03 

9 

3/6/86 

20:05 

59 

8.3 

3.91 

10 

3/6/86 

22:13 

33 

5.0 

9.42 

11 

8/29/86 

23:01 

11 

2.5 

24.33 

12 

9/23/86 

21:21.5 

14.5 

3.2 

17.03 

13 

10/18/86 

20:17.5 

16.5 

4.5 

13.47 

14 

11/4/86 

16:54 

13 

5.3 

17.64 

15 

11/15/86 

1:46 

16 

2.8 

14.29 

17 

11/24/86 

20:29 

8 

3.2 

34.61 

18 

12/21/86 

22:58 

12 

3.2 

30.70 

19 

12/23/86 

21:53 

18 

3.2 

19.43 

20 

1/18/86 

0:42 

16 

2.8 

23.02 

Event  date  and  start  time  were  based  on  GOES  5  observations, 
time  delay  is  measured  from  the  wave  onset  time  at  GOES  5  to  the 
wave  onset  time  at  GOES  6. 


